The Fibonacci quantum walk and its classical trace map 
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We study the quantum walk in momentum space using a coin arranged in quasi-periodic sequences 
following a Fibonacci prescription. We build for this system a classical map based on the trace of the 
evolution operator. The sub-ballistic behavior of this quantum walk is connected with the power-law 
decay of the time correlations of the trace map. 
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I. INTRODUCTION 

The quantum walk (QW) on the line is a subject that 
has drawn much attention in quantum computing [1]. It 
has been introduced by 0, Q j as counterpart of the clas- 
sical random walk. Due to the success of the classical 
random walks in the description of many classic models 
of the nature, the possibility that future quantum algo- 
rithms will be based on the QW has attracted the atten- 
tion of researchers from different fields. Additionally, as 
a result of recent advances in experimental techniques to 
manipulate quantum states has led the scientific commu- 
nity to think that quantum computers could be a reality 
in the near future. Anyway the study of the QW sub- 
jected to different sources of decoherence is a topic that 
allows to analyze and verify the principles of quantum 
theory and has been considered by several authors 

One of the most striking properties of the QW is its 
ability to spread over the line linearly in time as char- 
acterized by the standard deviation a(t) t, while 
its classical analog spreads out as the square root of 
time [cr(t) ~ t^/^]. Recently, in ref. Q the QW was 
generalized introducing two coin operators arranged in 
quasi-periodic sequences following a Fibonacci prescrip- 
tion. This "Fibonacci QW" leads to an unexpected 
sub-ballistic wave function spreading, as shown by the 
power-law tail of the standard deviation [cr(i) ~ with 
0.5 < c < 1]. More recently [1, 0], we have studied the 
QW subjected to noise with a Levy waiting-time distribu- 
tion [Sj], finding that the system has also a sub-ballistic 
wave function spreading instead of the known ballistic 
growth. The Fibonacci sequence gives rise to a rich be- 
havior also in other systems; it is the cause of sub-ballistic 
behavior both for the quantum kicked rotor in resonance 
and for tight-binding electrons [13, [U • On the other 
hand Sutherland p^ . twenty years ago, studied a very 
simple two dimensional system of a single particle with 
spin 1/2 subjected to quasi-periodic external field pulses, 
showing that the dynamics of the system is determined 



by a three dimensional "trace map" . In that work, he 
found three regions of the parameters space character- 
ized by qualitative differences of the decay of the time- 
correlation function: slower than a power law, as a power 
law, and faster than a power law. 

Here we study the Fibonacci QW through the dynam- 
ics of its trace map. We show that the sub-ballistic be- 
havior of the QW is connected with the power- law decay 
of the time-correlation function of the trace map. 



II. FIBONACCI QUANTUM WALK 

The dynamics of the Fibonacci QW will be generated 
by a large sequence of two time-step unitary operators 
Ui and U2 for each time t. Given the two initial values 
of the succession Ui and U2 the sequence is obtained, as 
was done in previous works [1, [H, [3 1 using the rule 



Uk 



+1 



UkUk- 



(1) 
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To obtain the operators Ui and U2 we develop in some de- 
tail the free QW model. The standard QW corresponds 
to a one-dimensional evolution of a quantum system (the 
walker) in a direction which depends on an additional 
degree of freedom, the chirality, with two possible states: 
"left" \L) or "right" \R). Let us consider that the walker 
can move freely over a series of interconnected sites la- 
beled by an index n. In the classical random walk, a 
coin flip randomly selects the direction of the motion; in 
the QW the direction of the motion is selected by the 
chirality. At each time step a rotation (or, more gener- 
ally, a unitary transformation) of the chirality takes place 
and the walker moves according to its final chirality state. 
The global Hilbert space of the system is the tensor prod- 
uct Hs IE) He where Hg is the Hilbert space associated to 
the motion on the line and He is the chirality Hilbert 
space. 

We shall consider only unitary transformations which 
can be expressed in terms of a single real angular param- 
eter 9 [3, [H, [i3|. Let us call T_ (T+) the operators 
that move the walker one site to the left (right) on the 
line in Hs and let \L){L\ and |-R)(-R| be the chirality pro- 
jector operators in He- Then we consider free evolution 
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transformations of the form [14| , 

U,{9,) = {r_ ® \L){L\ + r+ ® \R){R\} o {/ ® if (0,)} , 

(2) 

where K{6i) — az& is an unitary operator acting 

on i?c, CTy and tJz being the standard Pauh matrices, and 
J is the identity operator in Hg. The unitary operator 
Ui{9i) evolves the state by one time step, 

\^it+l))^UM)\'i'it))- (3) 
The wave vector is expressed as the spinor 



bn{t) 



(4) 



where we have associated the upper (lower) component 
to the left (right) chirality, the states \n) are eigenstates 
of the position operator corresponding to the site n on 
the line. The unitary evolution for |^(i)), corresponding 
to Eq. ([3]) can then be written as the map 

a„(i + 1) a„+i(t) cos 9 + bn+i{t) sin 9 , (5) 
bn{t + 1) = a„_i(t) sin 9 — bn-i{t) cos 9 . 

To build the operators Ui and U2 we substitute in the 
previous expression 9 by 9i and 02 respectively. As men- 
tioned in the introduction these types of maps produce a 
sub-ballistic wave function spreading. 



III. TRACE MAP 

With the operators defined in the previous section, 
we shall study the Fibonacci QW connecting it explic- 
itly with the work of Sutherland [l^l. To this end, 
we make the spatial Fourier transform of the ampli- 
tude (a„(t), 6„(t))^ multiplying both sides of Eq. ([5]) by 
gz(0-7r/2)n^ with (f> G [— 7r,7r], and summing in the integer 
index n: 



where 



Unit) 



and 



, , . , e *"^cos 6* e sm 

M{(p,9) =1 id, ■ a i<h 
^ ' * sm 9 cos 



(6) 

(7) 
(8) 

(9) 



Thus in Fourier space, the dynamics of the Fibonacci QW 
is determined by the 2x2 unitary matrix M{(j),9). We 
call Ml and M2 the matrix M{(j>, 9) evaluated in (0i, 6*1) 



and {(j)2,92) respectively. The matrix M{(j>,9) can be 
rewritten employing the Bloch-sphere representation as 

M(0, 6*) = cos6lsin(/) / + «M.(T, (10) 

where / is the identity matrix, u is given by 

u=(sin 6* cos 0, sin 6* sin (/), cos 6' cos 0) , (11) 

and cr = {<Jxi cTy, az)- The vector u can point in any di- 
rection but for our purposes only the relative orientation 
between ui and U2 is important. If we choose (jji = tt/2 
then 



Ml = cos 9i I + i sin 9i (t„ 



(12) 



The values of the parameters 9i, 92 and 4>2 determine 
the initial state for the dynamics of the system. In the 
Fourier space, the system can be interpreted as a particle 
with spin 1/2 subjected to two types of pulses. Mi and 
M2, coming from two external magnetic fields [13]. The 
dynamics of the spin is specified by the order of the in- 
teraction pulses. In our case the pulse Mi is apphed first, 
then the pulse M2 is applied and the following pulses are 
determined by the rule 



Ml 



k+i 



= MkMk-i. 



(13) 



The pulse Mi is a rotation with axis pointing in the pos- 
itive J direction, then it is completely determined spec- 
ifying 



xi = TrMi/2 = cos6li. 



(14) 



where TrA is the trace of the matrix A. The pulse M2 
can point in any direction but it is determined specifying 
02 and (j32', or equivalently by 

yi = rrM2/2 = sin02COs6l2, (15) 
zi = rr(Af2Mi)/2 = sin02Cos(02 + 0i). (16) 

Starting from the ordered pair (Mi,M2), called stan- 
dard alignment, any other orientation {Mk, Mk+i) can 
be thought as a rotation of SU{2). That is 

(Mfc, Mfe+i) = R (Ml, M2) R-\ (17) 

where R is the rotation operator. The vector rj, — 
(xfe, j/fc, Zk) is defined as 



1 



rk 



{TrMk, TrMk+i, TrM, 



k+2) 



(18) 



Now, the dynamics of the system can be visualized 
through the iteration of the trace map (xk, yk, ^fe) =^ 
(xfc+i, yk+i, Zk+i) with initial condition {xi, yi, zi). 
Using some properties of these matrices, Kohmoto et al. 
[3] showed that the vector satisfies the following non- 
linear map 



{xk+i, z/fe+i, Zk+i) = [yk, Zk, 2ykZk-Xk). 



(19) 
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FIG. 1: Poincare sections for the two hemispheres of the trace 
map. Arbitrary units are used. For C=-0.99 a) back and b) 
front hemispheres. For C=-0.7 c) back and d) front hemi- 
spheres. 

This map has an invariant given by 

C = xl+yl + zl-2xkykZk-l = - (sin 6li cos ^i2)^ (20) 

The value of the invariant is obtained by substitution of 
the initial condition in the generic expression of C. It 
also follows from Eq. ^ that -1 < C < 0. In Fig. [1] 
and Fig. [2] we show a sequence of Poincare sections of 
the trace map in the plane {x,z) for increasing values 
of C. This sequence shows that the size of the sections 
grows with the value of the invariant C. As the invariant 
surface has the topology of a sphere we represent sep- 
arately both hemispheres corresponding to positive and 
negative values of y. The typical orbits of the map are 
presented over both hemispheres. For each value of C the 
left (right) figure is the back (front) hemisphere. Both 
hemispheres look similar, as if there existed a symmetry 
operation between them. Each hemisphere presents sev- 
eral axes of symmetry and some elliptic and hyperbolic 
points. 

In Fig. [Tla,b for C — —0.99 there appear to be no 
chaotic regions for the scale shown. It is easy to identify 
four elliptic points and one hyperbolic point with differ- 
ent winding numbers. Therefore, for values of C near —1 
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FIG. 2: Poincare sections for the two hemispheres of the trace 
map. Arbitrary units are used. For C=-0.53 a) back and b) 
front hemispheres. For C=-0.5 c) back and d) front hemi- 
spheres. 



there appears to be a second integral of the motion, at 
least to the accuracy of these plots, and the system seems 
to be integrable. 

In Fig. [l]c,d for C — —0.7, the largest resonances are 
clearly seen, but some of the periodic orbits are pushed 
to the periphery. The neighborhood of the central hy- 
perbolic fixed point begins to be chaotic but the plot is 
still dominated by the KAM tori. 

In Fig. [21a,b for C = —0.53, the central hyperbolic 
fixed point has a chaotic sea but there are still some el- 
liptic points stretching diagonally. 

Fig. [He, d for C = —0.5, is similar to the previous case 
but now the each elliptic zone breaks into two pieces and 
the chaotic sea increases. 

It is interesting to study the trace map when ^2 — ''■/2, 
and then C — 0. As showed in Fig. [3] the second integral 
of motion appears to have been totally destroyed and 
the integrability of the system is lost. In this case the 
invariant surface has two accumulation points. In the 
other extreme, when C — —1, the solution of the trace 
map is only one point. 

For intermediate values of C chaos sets in and large 
regions contain a chaotic sea of trajectories. There is 
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no way to predict the future evolution of a trajectory in 
these regions, but it is possible to use the techniques of 
stochastic theory to determine the statistical behavior of 
trajectories. If no stable islands or cantori exist, the diffu- 
sion process is like that found in simple Brownian motion. 
However, if stable islands and cantori exist, the diffusion 
process becomes much complex. Here we see that there 
are many stable islands embedded in the chaotic region. 
The darker shading of trajectories along some lines indi- 
cate that these regions are bracketed by cantori. During 
the wandering in the chaotic sea the trajectory could get 
trapped in the neighborhood of a stable island for arbi- 
trary lengths of time. The trapping of trajectories by the 
stable islands gives rise to power-law decay of trajectory 
correlations rather than exponential decays as would be 
expected for a purely random diffusion process. Suther- 
land 12] has computed the correlation function for the 
trace map. Our equations are the same as in Ref. p^ . 




FIG. 3: Front hemisphere for the trace map whit C=0, the 
back hemisphere looks the same up to a rotation. 

but our initial conditions are chosen obeying the initial 
conditions of the Fibonacci QW. In both models: i)the 
identification of the most obvious elliptic and hyperbolic 
-ftT-cycles are similar, ii) he relation between the rotation 
of the standard alignment (Mi, Mi) and the ii'-cycles 
of the trace map are common, iii) the trace map dy- 
namics determines the rotation dynamics. Therefore we 
can translate the results of the Sutherland model to the 



Fourier transform of the Fibonacci QW. Then, here also 
three regions of the parameter space are characterized by 
the qualitative differences in the decay of the time corre- 
lations. These regions are: i) near a stable elliptic cycle 
where the decay is slower than a power of time, ii) near 
a hyperbolic cycle where the decay is a power of time, 
iii) in the chaotic region where the decay is faster than a 
power of time. 

IV. CONCLUSION 

In this work, we have established a trace map for the 
Fibonacci QW. This map has an integral of the motion 
C that determines a three dimensional surface where the 
trajectories lie. The trajectories depend on the param- 
eter C, when this parameter goes to —1 the system is 
integrable, but when it goes to the trajectories become 
chaotic. For intermediate values of C the trajectories 
have a rich behavior where it is possible to see on the 
(cc, z) plane elliptic orbits, hyperbolic orbits, KAM tori, 
cantorus, etc. It is in this zone where the autocorrelation 
function has a time power-law dependence. This devia- 
tion from exponential decay appears because the stable 
islands in the chaotic sea can trap trajectories for long 
intervals of time. This stickiness of the islands seems to 
be due to a network of cantori and higher order stable 
islands chains which surround the stable islands and ex- 
tend into the chaotic sea. It is possible to show, as in Ref. 
[TtI , that the sub-ballistic increase in the variance of the 
Fibonacci QW is a direct consequence of the coherence 
of the quantum evolution. Then, not all the trajecto- 
ries of the trace map contribute with the same weight in 
the Fourier antitransformation of Eqs. (|7l8p . The more 
regular behavior has more weight than the chaotic trajec- 
tories and this leads to the power-law decay. Therefore, 
the power-law behavior in the trace map affects the wave 
function in such a way to obtain the sub-ballistic behav- 
ior of the standard deviation. However, it is clear that 
more work will be necessary to understand these mecha- 
nism more deeply. 
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